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Several no-go theorems showed the incompatibility between the locality assumption and 
quantum correlations obtained from maximally entangled spin states. We analyze these 



no-go theorems in the framework of Bohm's interpretation. The mechanism by which non- 



local correlations appear during the results of measurements performed on distant parts 



of entangled systems is explicitly put into evidence in terms of Bohmian trajectories. It 
is shown that a GHZ like contradiction of the type+l=-l occurs for well-chosen initial 
q ■ positions of the Bohmian trajectories and that it is this essential non-classical feature that 



makes it possible to violate the locality condition, 
keywords: Bohmian dynamics, non-locality 

Introduction. 



Several no-go theorems show that it is impossible to conciliate locality and quantum me- 
chanics. The essence of these theorems is that it is impossible to simulate the results 
of observations carried out on distant parts of an entangled quantum system in terms 
of a local common cause mechanism. This impossibility can be expressed through the 
violation of an inequality f|, [15|] or, even "stronger", of an equality |10|, [HJ. Besides 
it is well known that the Bohmian interpretation according to which particles have at 
any time a well-defined position || makes it possible to simulate all the quantum pre- 
dictions for massive particles in the non-relativistic regime. Therefore, it is interesting 
to understand better how non-locality appears in the Bohmian picture. It is well-known 
that for entangled systems the velocity associated to the Bohmian trajectories can be 
non- locally influenced []5], |J. The goal of our paper is to analyse in detail the nature 
of the Bohmian trajectories in the specific situations encountered in the formulation of 
three well-known no-go theorems, Bell's theorem Jj]], Mermin's theorem |13| and GHZ's 



theorem [1C], in which respectively two, three and four spin 1/2 particles are assumed to 
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be prepared in maximally entangled states. The knowledge of the initial positions of these 
particles is sufficient in Bohm's formulation in order to predict the whole trajectories and 
also the results of local spin measurements performed with Stern-Gerlach devices. The 
specification of the initial positions can thus be considered as a common cause mechanism 
which makes it possible to explain all the quantum correlations. In accordance with the 
no-go theorems and with the properties of the Bohmian dynamics, this common cause 
mechanism may however not be considered to be local. We shall show explicitly in the 
first section of this paper how, for well chosen initial positions, a GHZ like contradiction 
occurs for a situation during which Bell's inequalities get violated. This result confirms 
Hardy's formulation of Bell's inequalities [[II] according to which the violation of Bell's 
inequalities is proportional to the probability that a GHZ like contradiction of the type 
+ 1=-1 occurs flQfl . We shall also show explicitly in the second section of this paper how 
in the situation considered in the no-go theorems of Mermin [T3] and Greenberger, Home 
and Zeilinger |10[ the GHZ like contradiction occurs for all the initial positions. 



1 Bell's theorem in a Bohmian description. 

1.1 Simultaneous Stern-Gerlach measurements performed on a 
pair of entangled particles. 

In the situation described in Bell's theorem ffl], two particles are prepared in an entangled 
spin state and spin measurements are performed simultaneously on each particle. The 
full wave-function associated to such a state can be put in the form: 

\V(r L , r R , £)) = 

a^J ++ {r L , r R , t) ■ | +) L ® | + } R + H+-{r L , r R , t) ■ | +) L ® | - ) R 

+ c^_+(r L , r R , *) ' I - >l ® I + >r + #— (rx, r R , t) ■ | - ) L ® | - ) R (1) 

where a, b, c, d are complex constants representative of the spin entanglement, the indices 
L and R represent arbitrary spatial reference frames in two distant regions, the Left and 
Right regions, r L / R represent the position vector expressed in these frames, and | +/ — ) L / R 
represent up/down spin states along the (not necessary parallel) Z-axes of these frames. 
We shall assume for convenience that Stern-Gerlach devices are disposed at the same 
distance from the source, along the axes of propagation Xl/r of the particles, that the 
point of penetration inside them coincides with the origins of the spatial frames, and with 
the origin of time, and that the initial state consists of a pair of Gaussian shaped particles 
of mass m propagating at a speed vq : 

* (rL ' P *' t]) = (2 V ^L ) 6 / 2 ex P C r l~2 R )exp{ik {x L + x R )). (2) 

When the Zl/ r axis is parallel to the L/R magnetic field, it can be shown that the 
equation of evolution is separable into the coordinates (xl, Vl, Zl, xr, Ur, z r ). For reasons 
of convenience, we shall also assume that the Stern-Gerlach set-ups are similar in the Left 
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and Right regions. Then, it is easy to show that the solution of the equation of evolution 
given in the appendix 2 can be expressed in terms of the solution of the single particle 

case: 

\V(r L , r R , t)) = 

aip+ lOl, t) ■ ip+ R (r R , t) ■ | + ) L <g> | + ) R + bvp + L (r L , t) ■ ^- R (r R , t) ■ \ + ) L ® | - ) R 
+ a/)- L (r L , t)-ip + R (r R , t)-|-) L ®| + ) R + #- l(t l , t)-if>- R (r R , t) ■ | -) L <g> | -) R , (3) 
with: 

'0+/- l/r(*l/r, t) = l/r (x L /r) • l /r(Vl/r, t) ■ l/r( z l/r) (4) 

where 4>+/-i an d V'+y- are ^ ne solutions of the single particle case that were origi- 

nally derived by Bohm 0. This derivation is reproduced in detail in the appendix 1. 



1.2 Bohm's interpretation. 

Bohm's interpretation is, in summary, the following: whenever we can associate to the 
equation of evolution of a quantum system a conservation equation of the form dtp = 
div(J'), where p is a positive definite density of probability, and div(J) is the divergence of 
a current vector, we can interpret this density as a distribution of localised material points 
moving with the velocity — . It is thus possible to formulate a hidden variable theory for 
the system: it would consist of a spatial distribution of material points which initially 
coincides with the quantum distribution (given by p); these points move with a velocity 
equal to J j p. In virtue of the conservation equation, the spatial distribution deduced 
from this evolution coincides then for all times with the quantum distribution. According 
to de Broglie, all the measurements being, in last resort, position measurements, this 
hidden variable theory is, for what concerns practical purposes, equivalent to orthodox 
quantum mechanics. 

In the appendix 2, we show how to deduce such a conservation equation for the pair 
of particles considered here. From this equation of conservation, it is straightforward to 
deduce the conserved density (eqn 50): 



p{r L , r R , t) = (tf(r L , r R , t)|tf(r L , r R , t)) = 

\aip+ L (r L , t) -ip + R (r R , t)\ 2 + \bip+ L (v L , t) • V- r(t r , t)\ 2 + 

\af>- L (r L , t) ■ V>+ r (t r , t)\ 2 + |#- l{t l , t) ■ ^- R (r R , t)\ 2 (5) 

This density depends in general on the two positions (r^, r^). The 6-dimensional velocity 
(eqn 51) associated to the conserved density through the relation J = p- v can be splitted 
into two local 3-dimensional velocities vl, v r defined as follows: 

v L (r L , r R , t) = — (3f(V»+ L (r L , t)* V L ^+ L (r L , t))-(\a^ + R (r R , t)\ 2 + \b^ R (r R , t)\ 2 ) 
mp 

+ 3(^_ L (r L , t)* V L ^_ L (r L , t)) ■ (|c0+ R (r R , t)\ 2 + |#_ R (r R , t)\ 2 )) (6) 
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vr(*l, r R , t) = — (9f(^+ fl(r fl , *)* Vfl^+ «( r ii, *))-(l a ^+ l^l, ^)| 2 + |c^- t)| 2 ) 
mp 

+ 3(^_ fl (r fl , *)* V fi ^_ R (r R , *)) • (\bil> + L (r L , t)\ 2 + |#_ L (r L , t)| 2 ) (7) 

where Vi(m is the gradient on the spatial Left (Right) coordinates, and $>(z) the imagi- 
nary part of z. In general, these velocities contain a non-vanishing contribution depending 
on the position of the particle situated at the other side of the source, and are thus influ- 
enced by what happens in a region from which they are separated by a space-like distance. 
In the formulation of Bell's theorem, the particles are assumed to be prepared in the so- 
called singlet state for which a = d = ■ sm(|), and b = — c = ■ cos(|) where 9 is 
the angle between the magnetic fields in the Left and Right Stern-Gerlach devices. When 
these axes are parallel, we recover the standard expression: 

|*(r L , r R , t)) = 

^ + L (r L , t)-i,_ R (r R , t )-\+) L ®\-) R - -J=V- l{t l , t)-i> + R (r R , t)-|-) L ®| + ) R (8) 

It is worth noting that when the axes of the devices are parallel, the results of spin mea- 
surements are always opposite. This perfect correlation is due to the fact that the singlet 
state is maximally entangled, and that its expression is invariant under a simultaneous 
rotation of the axes of quantization, when both axes are parallel. The specification of a, b, 
c and d allows us to express explicitly the Bohmian dynamics during the passage through 
the Stern-Gerlach devices. We obtain, after some lenghty but simple computations, the 
following velocities for the Left as well as for the Right x and y components: 

dec Jv^ 

t - rrw* (10) 

(where k = 2 m5 2 r ex P resses the diffusion inherent to the Schrodinger equation). The 
solutions of these dynamical equations are: 



x 



v t + x Vl + k 2 t 2 (11) 



y = y Vl + k 2 t 2 (12) 

In the appendix 1 we show that, during a typical Stern-Gerlach experiment, the diffusion 
is negligible during the time of flight through the Stern-Gerlach. The trajectories along 
these axes are thus essentially free translations, and present no special interest for the 
understanding of the measurement process. 

In the singlet state, the dynamics of the ^-components is given by: 

dz L k 2 

~dT = i + kH 2 tZL + 

I fit 2 (z L + z R ) s , ()t 2 (z L -z fl K , fit 2 (~ z L + z r) \ ( fit 2 ' i- z L ~ z r) \ 

_l_ g e 1 + k2t ' 2 2 + c e 1 + k2t ' 2 2 — c e 1 + k ' 2t ' 2 2 — s e 1 + k ' 2t ' 2 2 



, fit 2 ( z l+ z r) \ I lit 2 ( Z L - Z R.) \ ( fit 2 (- Z L + Z L) \ I fit 2 

_|_ S 2gV 1 + k 2 t 2 2 > _|_ C 2gV! + fc 2j2 2 > _|_ C 2g4 + fc 2 t 2 2 > _|_ g2 g 1 . 1 + fe 2 t 



fit 2 (-zr.~ z l) \ 
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I fit 2 (z L + z R ) , , pt 2 (z L -zrK / fit 2 (~ z l + z l) \ I 0t 2 {-*R-*l) \ 

2 p y l + k^i 1 2 > _ p^p^ l + fc2 t 2 2 ^ I ptp^i + fc^ t 2 2 J _ + fc^ t 2 2 ^ 



•{ rT ^«t 3 + 2at} (13) 
k 2 

~dT = i + k 2 t 2 tZR 4 

+ s z e'' 1 + fc2t2 2 ' — c z e Kl + fc2 ' 2 2 + c e 1 + fc2 * 2 2 ' — s z e Vl + fe2i 

/ lit 2 ( z L+ z r) \ I lit 2 i z L - z r) \ I fit 2 (- Z L + Z R) \ I lit 2 (_Z±L ~ Z R.) \ 

-\- s 2 e 1 + fc2j2 2 + c 2 e 1 + fc2 * 2 2 + AhW 2 J + s 2 e 1 + fc2j2 2 

■ { rrw at3 + 2at} (14) 

where a = ai/x/2m, /3 = a\^ijmbr\, c = cos(|) and s = sm(|). 

Numerical computations (confirmed by a careful evaluation of the terms present in 
the last equations) show that, when the two magnets act simultaneously, there exist 
four attractors for the trajectories, corresponding to the results of spin measurements 
(upL, upn), (upL, downn), (downi, upn) and (dowrii, dowriR). During the passage 
through the device, the trajectories can be shown to fall very quickly in the basin of one 
of these attractors. For instance, in the attractor basin of the (upi, upr) outcome, the 
( fJt2 ., < - ZL + Zr) ) 

exponential factor e 1 + fc2 ' 2 2 will dominate all the other exponential factors, and 
the z-velocities will obey the following equations: 



dz R k 2 — k 2 

tz R +1 



dt i + 1 ' { TTT^ at +2at} 

In virtue of the fact that, in a well-conceived device, we can neglect the diffusion terms 
relatively to the classical drift terms, we can take k to be equal to in the previous 
equation, which gives the classical expression for the velocities (2a is the magnitude of 
the classical acceleration undergone inside the magnetic field of the device): 

dz r . 

= + 2at 



dt 

dz R 



= +2at 



dt 

Similarly, in say the attractor basin of the (down L , upr) outcome, we shall get: 

dzj, 

= -2at 
= +2at 



dt 

dz R 



dt 

This means that after a short time, the trajectories become and remain quasi-classical. 
The choice of the attractor basin is deterministically speficied by the initial positions. 
It is worth noting that when the initial positions lie in between the regions that are 
characterised by this quasi-classical behavior the deflection will be quite less important. 
These positions can be associated to the tails of the outgoing wave-packets and their weight 
is negligible, because the diffusion of the Gaussian wave-packets is negligible relatively to 
the average deflection undergone during the passage through the Stern-Gerlach device. 



5 



1.3 GHZ like contradiction in the two particles case. 

In last resort, the final outcome is thus "nearly always" unambiguously determined by the 
initial positions zl and zr. Let us assume that we choose to change the direction of the 
magnet in the Right device. Then, the final outcome is determined by the initial positions 
zl and z' R . z' R depends on yn and zr and on the angle between the old quantization axis 
Zr and the new one Z' R . It could happen that for the same initial position coordinates 
Z L-, Vri z Ri the attractor basin of the Left trajectory changes due to the change of zr into 
z' R . This means that the choice of the basis of detection in the Right region can non- 
locally influence the outcome of the measurement in the Left region. Therefore although 
the results of all possible measurements can be deterministically foreseen on the basis 
of the knowledge of the initial positions, which constitutes a common-cause mechanism, 
this mechanism is non-local. This helps to understand how Bohmian mechanics makes 
it possible to violate Bell's inequalities. The probability of an outcome can be obtained, 
after integration of the Bohmian velocities, by performing a weighted average on the 
predictions associated to this outcome for particular initial positions. The weight is equal 
to the initial probability of presence, a Gaussian distribution in our case. In accordance 
with Bohm's interpretation, the probability obtained so is exactly the same as the standard 
quantum probability. In [0] , the attractor basins of the outcomes related to different non- 
compatible spin measurements in the Left and Right region were explicitly determined 
for a well-chosen initial position and the following result was shown. If we choose four 
directions {Zl, Z' L) Zr, Z' r ) such that the directions Zl, Z' l , Z' r are coplanar and are 
all separated by angles of 120 degrees, while Z' L and Z R are parallel, that the system of 
two particles is prepared in the singlet state, and that the initial coordinates of position 
of the pair are defined as follows: (y L , z L )(t = 0) = 10" 3 cm ■ (y/3/2, 1/2), and 
{Vri z r){P — 0) = 1.1 • 10~ 3 cm • ( Vo/2, 1/2) then we shall observe, after integration 
of the Bohmian velocities, the following outcomes: (upz L , upz R ), (upz' L , Wz' R ), and 
(down Z ' L , upz R )- R the locality assumption was valid, the results of the measurements 
ought to be predetermined, before the measurements take place. This can be shown by 
an EPR like reasoning |§ which is, roughly summarised, the following: -in the singlet 
state, the value of the spin component along an arbitrary axis of reference in the Left 
region can be deduced from the observation of the corresponding value in the Right 
region because the singlet state exhibits perfect (100 %) correlations; -if we assume that 
distant measurements may not influence the local spin values, then these values existed 
necessarily before their measurement took place. This means that the local value (in 
the Left region but also in the Right region by a similar reasoning) of the spin along an 
arbitrary direction ought to be predetermined before the measurement. 

Obviously, if the locality assumption was valid, we could deduce from the predictions 
relative to the two first joint measurements ((upz L , upz R ), (upz' L , upz')), that the third 
result ought to be (upz> L , upz R ) too, which is not the case here (we find (dowriz> L , upz R ))- 
In [|7| we showed how Bohmian non-locality is related to the non-factorisability condition 
postulated by Bell as a consequence of the locality condition. For the purposes of the 
present paper, we went further. We also integrated the velocities when the magnets are 
parallel to the direction (Z L , Z' R ), and we found the outcome {upz L , upz> )• By assigning 
a value +1 to the outcome spin up and -1 to the outcome spin down, we arrive to a GHZ 
like contradiction of the type +1=-1 that we shall study in detail in the next section. 
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This can be done by considering the product of the outcomes of Z L fl Z R and Z' L fl Z' R 
at one side and the product of the outcomes of Z' L fl Z R and Z L fl Z' R at the other side. 
If the outcomes associated to the four experiments were predetermined, these products 
ought clearly to be equal but for the choice of initial position considered here they are 
respectively equal to +1 and -1, which can be put in the paradoxical form +1=-1. 

It can be shown j7j that when the four directions (Zl, Z' l , Z r , Z' r ) are coplanar 
and are all separated by angles of 120 degrees, while Z' L and Z R are parallel, and that 
the system of two particles is prepared in the singlet state, Bell's and Clauser-Horne's 
inequalities are violated. It can also be shown that a GHZ like contradiction occurs only 
for a fraction of the initial positions that differs significatively from unity. For instance 
when the initial coordinates of position of the pair are defined as follows: (i/l, ZL){t = 
0) = I0~ 3 cm ■ (V3/2, 1/2), and (to, z R )(t = 0) = -1 ■ l(T 3 cm ■ ( y/3/2, 1/2) 
then we shall observe, after integration of the Bohmian velocities, the following outcomes: 
(upz L , downz R ), (upz' L , down Z ' R ), (upz L , down Z ' R ) and (upz' L , down ZR ). In this case, no 
GHZ contradiction occurs, and the results of local measurements do not depend on which 
measurements are performed in a distant region. These results confirm Hardy's analysis 
who showed [TTI] that the degree of violation of Bell's and Clauser-Horne's inequalities is 



proportional to the probability of encountering GHZ like paradoxical situations. 



Besides, Peres [0 showed that for two particles in the singlet state such a paradoxical 
situation also characterised the quantum average values of well-chosen observables. The 
interpretation of Peres's paradox in the framework of Bohm's theory was already studied 
by Dewdney in || so that we shall not repeat his work in the present paper. Furthermore, 
it is impossible to express Peres's paradox in terms of local Stern-Gerlach measurements 
only because it involves non-local measurements. Therefore, presently, there is no explicit 
formulation of this paradox in terms of Bohmian trajectories. Nevertheless, such a formu- 
lation is possible in principle || . It is worth noting that the implications of Bohm's theory 
regarding non-locality were already recognised by Bohm and Bell, and that the elucida- 
tion of certain no-go theorems was qualitatively developed in the past (see for instance 
the references || U and references therein). Computer simulations of the dynamics inside 
the Stern-Gerlach devices (with two particles and parallel magnets only) can be found 
in |J, as well as a qualitative discussion of their connection with Peres's paradox. The 
analytical expression of Bohm's velocities in the general two particle case (with parallel 
and non-parallel magnets) was given in as well as the connection between Bohmian 
non-locality and non-factorisability. Nevertheless, the occurence of a GHZ like contradic- 
tion in the two particles well as the explicit study of the topology of the attractor 
basins for the three and four particle cases and their connection to GHZ like contradic- 
tions were not, to the knowledge of the authors, published elsewhere. This brings us to 
the next chapter. 



2 Mermin and GHZ no-go theorems in a Bohmian 
description. 

The two paradoxes that we shall study in this section are of the following form: for well- 
chosen local observables and quantum states, the quantum correlations are such that if 
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the results of the local measurements are fixed by advance and do not depend on the 
measurements that are performed in distant regions we come to a contradiction of the 
type +1=-1. We showed in the previous section that in the two particles case such a 
contradiction could occur for well-chosen initial positions. However, after averaging over 
the initial positions of the particles, for which the contradiction occurs sometimes but 
not always, the paradox is expressed by the violation of an inequality in accordance 
with Hardy's analysis |fTI |. The novelty of Mermin's and GHZ paradoxes |10|, |13[ is that 
the contradiction does not take the form of the violation of an inequality but well of an 
equality. As we shall now show, the meaning of this result, interpreted in terms of Bohm's 
trajectories is that all the initial positions lead to the same paradoxical situation, and not 
only some of them as in the two particles case. 



2.1 Mermin's no-go theorem: GHZ like contradiction in the 
three particles case. 



In Mermin's theorem [fl3 1 , three particles are prepared in a maximally entangled spin 
state of the form: 

|*(ri, r 2 , r 3 , t)) = 
i( r i» *) " ^+ a(r 2 , t) ■ 3 (r 3 , t) ■ | + ) x <g> | + ) 2 <g> | + ) 3 

- t) ■ V>- a(r 2 , t) ■ V>_ s(r 3 , t) ■ \ - ) x <g> | - ) 2 <g> | - ) 3 , (15) 

The particles are sent along three coplanar directions separated by an angle of 120 degrees 
and | +}; (| — }.) represents a spin up (down) state along the direction of propagation of the 
zth particle. Four "product" observables each of which is associated to a simultaneous spin 
measurement performed on the three particles are considered: E* • Y? y ■ Y? y , • Y? x ■ E^, E* ■ 
Y? y -Y? x and E^-E^-E^, where EJ. represents a Stern- Gerlach mesurement performed on the 
zth particle with the magnet parallel to the X direction which is orthogonal to the plane 
of propagation of the particles and E^ represents a Stern-Gerlach mesurement performed 
on the ith particle with the magnet parallel to the Y$ direction which is orthogonal to 
the direction of propagation of the ith particle and to X. The E matrices possess two 
eigenvalues, +1 and -1. The eigenvalue +1 corresponds to the outcome spin up and -1 to 
the outcome spin down. |^(ri, r 2 , r 3 , t)) is simultaneous eigenstate of the four product 
observables for the eigenvalues +1, +1, +1, —1. This imposes severe constraints on the 
results of local spin measurements. For instance, without performing the measurement of 
E^., we can deduce what would be the outcome observed during this measurement from 
the results of the measurements of Y? x and E^. By an EPR like reasoning similar to the 
one described in the previous section, the locality assumption implies that the outcome of 
the measurement of must be determined before the measurement took place and does 
not depend on which measurement is performed in distant regions on the other particles. 
A similar reasoning holds for the measurement of E^, E^, E^, E^ and E^. Let us denote 
S l x and S y the prediction assigned to the measurement of H l x and E^; the prediction +1 
corresponds to an "up" deflection inside the Stern-Gerlach device, and -1 to a "down" 
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deflection. The constraints imposed by the properties of the maximally entangled state 
in which the particles are prepared can be expressed as follows: 

sl-s 2 y -s 3 y = +1 

Sy • S% • Sy = +l 

Si -Si -Si = -1. (16) 
The product of the three first equations gives: 



(^) 2 • (S 2 y Y • (^ 3 ) 2 • Sl ■ S 2 X ■ S 3 X = +1 (17) 
S l y being equal to + or - 1, we get that 

Si- Si- Si = +1, (18) 

which together with the fourth equation leads to a contradiction of the type We 
encountered already such a contradiction in the previous section, and showed that it could 
be explained in terms of the non-local properties of Bohmian trajectories. This is also 
true in the present case. By a straightforward generalisation of the two particles case, it 
can be shown that the Bohmian velocities during the measurement of one spin component 
along X, one component along Y and another one along Y' are the following: 

dx k 2 

tx + 



dt 1 + k 2 t 2 

* 2 {x + y + y') \ , fit 2 (x-y-y') ^ , fit 2 i-x-y + y') ^ , fit 2 (-x + y-y 1 ) 



i pt \x t- y -i- y j \ / pi yj= — y — y) \ i pi \-x — y -r y >_ \ i pi 

-f e l i + k 2 t 2 2 > _|_ ^rrw 2 > _ e <-i + fe 2 t 2 2 > _ e \\ + k 2 t -z ■ 

f fit 2 (x + y+y 1 ) , , fit 2 (x-y-yl) s , fit 2 (-x-y + y 1 ) -, , fit 2 (-x + y-y') -, 

-\-e 1 + k ' 2 *' 2 2 + + e 1 + fe2 ' 2 2 e 1 + fc2 * 2 2 + e 1 + fc2i2 2 

■i 1 ' ^2 at " + ^ ( 19 ) 

dy k 2 



dt 1 + k 2 t 2 



ty + 



, (it 2 (x + y+y') s , (it 2 (x-y-y') , , fit 2 j-x - y + y 1 ') ) , fit 2 (-x + y-y') , 

_l_ e 1 + fe2 * 2 2 — e 1 + fc2 * 2 2 — e 1 + fc2t2 2 + e 1 + fc2j2 2 

/ /?t 2 (z + y + j/K / gt 2 (a: - v - »') \ / fit 2 (-x-y + y') , , fit 2 (-x + y-y 1 ) , 

-|- e 1 + fc2 * 2 2 + e 1 + fc2 * 2 2 + e 1 + fe2 * 2 2 + e 1 + fc2 * 2 2 

-A; 2 



^T~Tk^ at3 + 2at} (20) 

dy' A; 2 



dt 1 + Pt 2 

(3: - y - y') , 1 

-)-g v l + k 2 t 2 2 I — e V 1 + fc2 * 2 ' 2 + e 1 + fe2 * 2 2 ' — g v l + fe 2 t 

/ fit 2 (x + y + y') , , fit 2 (x-y-y') , . /it 2 (-x-y + y') , , fit 2 j-x + y - y') , 

-)- e l l + fc2 t 2 2 J _|_ gVi + fc 2 t 2 2 ' _|_ g 1 -! + fc2 t 2 2 -> _|_ g^l + fc2 t 2 2 J 



, (it 2 (x + y+y') , , (it 2 (x-y-y') , , fit 2 (-x-y + y') , , fit 2 (-x + y-y') , 

V l + fe2 t 2 2 > _ p 4 + fc 2 t 2 2 > I (="-1 + fe2 t 2 2 - 1 _ p y- L -+~~k T t I 2 J 



{ rT ^«t 3 + 2at} (21) 
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They possess four attractor basins that correspond to the outcomes (+ + +), (H ) 

( h) and (- + -). During the measurement of the fourth observable • • E^, the 

Bohmian trajectories obey: 



dxi k 2 

"dT = i + m 2 tXl 4 

+ e' 1 + fc2 * 2 2 ' + e v 



/ /3t 2 (-zi -£2 ~ ^3) -l ( P t2 (-^1 +^2+^3) 1 f /3t 2 (+^1 — + ^3) 'i / fit 2 j+Xi + x 2 - x 3 ) ^ 

— p ( l + k 2 t 2 2 > _ p( 1 + k 2 t 2 2 > I g^ + fc2 t 2 2 J I p\\ + k 2 t 2 2 > 



I fit 2 (-x 1 -x 2 ~ x a) ) I fit 2 (-£2 + x 2 + 3:3) \ / /3t 2 (+11 - 3=2 + gaj 1 I 111 . ^ +3: 1 +a: 2 ~ ^3) ■) 

4. g"> 1 + fc 2 t 2 2 ^ _|_ g V l + fc2 t 2 2 J _|_ gU + fc 2 t 2 2 - 1 _|_ g"" 1 + fe2 t 2 2 



■{r^^ + M ( 22 ) 

^2 k 2 

tXo + 



dt 1 + k 2 t 2 

+ x 2 + x 3 ) , 

4- g v l + ft 2 ! 2 " 2 ' — g v l + k 2 t2 2 1 _|_ gV 



/ gg* (-3:1 -3:2 - x 3 ) . , /3t 2 (-3:1 + 3:2+3:3) 1 f g£L (+^1 -3:2+3:3) ^ / fit- 1 (+gi +3: 2 -3:3) , 

3 l l + fc 2 t 2 2 J I + fc2 t 2 2 - 1 _ pVj T fc 2 t 2 2 - 1 I pVj T fc 2 t 2 2 - 1 



I 13 - ^ 1 / /3t 2 (-3;! + 3; 2 + 3:3) , , fit 2 ( + x± -x 2 + 313) \ , flt 2 ( + x 1 +x 2 - x 3 ) s 

-fe l i + k 2 t 2 2 1 _|_ g4 + fc 2 t 2 2 ' _|_ g^i + k 2 t 2 2 > _|_ g'-i + k 2 t 2 2 ) 

■{^« i3 + M (23) 
dx% k 2 

f ^ ~ 3:2 - ^ ~) I (it 2 (~x x +x 2 + x 3 ) , , fit 2 ( + 3?! - 3:2 + 3:3) \ / /it 2 ( + 3:! + 3: 2 - x 3 ) s 

— g l l + fc 2 t 2 2 -)- gM + k 2 t 2 2 J -)- g V l + k 2 t 2 2 > _ g^l + k 2 t 2 2 > 

, fit 2 (-x 1 -x 2 - x 3 ) - , fit 2 (.-xx +x 2 + x 3 ) . , flt 2 - 3:2 + 3:3) \ , flt 2 +3: 2 - 3:3) . 

-)- g l l + fc 2 t 2 2 ; 4- g V l + k 2 t 2 2 I 4- g^l + fc 2 t 2 2 J 4- g l l + k' 2 t 2 2 > 

■i x ~f kH , Q* 3 + 2at} (24) 

They possess four attractor basins that correspond to the outcomes ( ),(- + +),(+- 

+), and (+ + -)• Very quickly, the Bohmian trajectory will fall into one of these basins 
and follow a nearly-classical dynamics. For each choice of initial position (excepted for the 
initial positions that belong to the regions which separate two basins, which are regions of 
negligible weight that correspond to the tails of the Gaussian packets), the product of the 
predicted outcomes associated to the observables S^, Y? x and is equal to +1 when we 
consider their realisation during the three first experiments and to -1 when we consider 
the last one during which they are simultaneously measured. 



2.2 Greenberger, Horne and Zeilinger's no-go theorem: GHZ 
like contradiction in the four particles case. 

For what concerns the GHZ paradox [JTTJ, we shall firstly reformulate it in a simplified form 
that is closer to Mermin's formulation. Four particles are now prepared in a maximally 
entangled spin state of the form: 

|*(ri, r 2 , r 3 , r 4 , t)) = 
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-^=ip+ i(ri, t) ■ ip+ 2 (r 2 , t) • ^- 3 (r 3 , t) ■ ip- 4 (r 4 , *) • | + ) 1 <8> | + > 2 <8> | - ) 3 <g> | - ) 4 

- i(ri, *)-^_ 2 (r 2 , s(r 3 , *)-^+ 4(r 4 , *) • I - )i ® I - ) 2 ® I + ) 3 ® I + > 4 > ( 25 ) 

The particles are sent along four different coplanar directions and | + ) ; (| — ) i ) represents a 
spin up (down) state along the direction of propagation of the ith particle. Four "product" 
observables, each of which is associated to simultaneous spin measurements performed on 
the four particles, are considered: E* • T? x ■ Y? x ■ E^, ■ Y? x ■ E^ • H x , E* ■ Y? x ■ Y? x ■ E^ and 
E* ■ Y? x ■ Y? y ■ T,y, where EJ. represents a Stern-Gerlach mesurement performed on the ith 
particle with the magnet parallel to the X direction which is orthogonal to the plane of 
propagation of the particles and E* represents a Stern-Gerlach mesurement performed on 
the ith particle with the magnet parallel to the Y direction which is orthogonal to the 
direction of propagation of the ith particle and to X. |^(ri, r 2 , r 3 , r 4 , t)) is simultaneous 
eigenstate of the four product observables for the eigenvalues — 1, — 1, —1, +1. As before, 
this imposes severe constraints on the results of local spin measurements. By performing 
an EPR like reasoning similar to the previous ones, the assumption of locality implies 
that these results must be determined before the measurements take place and may not 
depend on which measurement is performed in distant regions on the other particles. 
Let us denote S x and S y the prediction assigned to the measurement of H l x and E* ; the 
prediction +1 corresponds to an up deflection inside the Stern-Gerlach device, and -1 to 
a down deflection. These constraints are expressed as follows: 

ql Q 2 q3 q4 1 

°i ' D i ' D x ' J i — 1 

ql Q 2 q3 q4 1 

' D x ' °y ' °x — 1 

ql Q 2 q3 q4 1 

°y ' °x ' °x ' °y — 1 

Si ■ S 2 X ■ Si ■ S A y = +1. (26) 
The product of the three last equations gives: 



Si ■ Si ■ Si ■ Si = +1, (27) 

which together with the first equation leads to a contradiction of the type +1=-1. The 
analysis of Bohm's trajectories is similar to the two and three particle cases so that we 
shall not repeat it entirely here. For instance, during the measurement of E* • Y? x ■ Y? x ■ E^ 
we get that the velocity of the first particle obeys the following equation: 

dx\ k 2 

~dt = 1 + kH* tXl + 

{-sh( - ^ k2t2 (-xi -x 2 - x 3 + x 4 )) - ^ k2t2 (-xi + x 2 + x 3 + x 4 )) 

fit 2 (it 2 

+ Sh( - k 2 f 2 ^ +Xl - ^2 + X 3 + X 4 )) + Sh( - k 2 f 2 ^ +Xl + %2 ~ %3 + X4))} / 

{ + ch( - ^ k2f2 (-xi -x 2 - x 3 + x 4 )) + ch( - ^ k2f2 (-xi + x 2 + x 3 + x 4 )) 
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+ ch (-, — !~TTZ? (+ X i - x 2 + x 3 + Xi)) + ch(- — 9 (+Xi + x 2 - x 3 + x 4 ))} 
1 + Art 2 1 + Art 2 

•{ rT ^at 3 + 2at}. (28) 

Similar equations are associated to the other velocities and to the other measurements. 
Only the factor that contains exponential terms changes from equation to equation. It 
is easy to guess the form of this factor because it obeys the following simple rules. To 
each outcome corresponds at the numerator and at the denominator as well a product of 
exponential factors of the form exp(( 1 _^ f fc2f 2 ( + /~ Xl+ /~ x 2+/-^3+/-x 4 ) ^ w j iere x . represents 
the projection on the axis parallel to the local ith magnet. The sign of Xj in the exponent 
represents the spin value asymptotically reached during the local measurement. The 
weight of this product is equal at the numerator and at the denominator as well to the 
probability to find the corresponding outcome. When the sign of Xi in the exponent 
is negative, the whole product will also be negatively weighted at the numerator when 
we consider the equation of evolution of Xj. Otherwise, and also at the denominator, 
the products will have a positive weight. This can be checked directly in the previous 
equation and is also true for the two and three particles cases. 

As in the three particles case treated in the previous section, the integration of the 
velocities leads to a paradoxical situation for any value of the initial position (excepted 
for some positions that belong to a region of negligible weight). 



3 Conclusions. 



In conclusion, we showed in this paper how the elucidation of several paradoxical situa- 
tions related to the question of non-locality can be realised in the framework of Bohm's 
interpretation and is finally reduced to the study of the topology of attractor basins of the 
Bohmian dynamics. One could object that this explanation is obvious being given that 
no satisfying relativistically covariant formulation of the Bohm's interpretation exists for 
situations in which more than one (entangled) particles are involved. This criticism is 
valid: in the present approach, we systematically considered the Schrodinger equation, in 
which a unique time appears, similar to the Newtonian, absolute, time. Now, a theorem 
due to Hardy Jl2| shows the impossibility to build a Lorentz-Invariant realistic theory 
that would mimick the predictions of quantum mechanics. In order to preserve realism 
it is therefore necessary to reintroduce a kind of absolute or etheric time, attached to 
an absolute or special frame of reference, which answers to the previous criticism. The 
interesting question is then, could the Bohmian interpretation help us to conceive an ex- 
periment that would reveal the existence of this quantum ether. The original Bohmian 
interpretation is too ad hoc to allow for such a possibility, but slightly modified versions of 
the interpretation (in which the randomisation of the hidden variable, here the position of 
the particle, according to the ip 2 distribution is not guaranteed for all times) would allow 
us to conceive such experiments |§. A realistic interpretation a la Bohm would then open 
the door to the conception of a quantum version of the Michelson-Morley experiment. 
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4 Appendix: Bohm's solution for the passage of one 
spin 1/2 particle through a Stern-Gerlach device. 



Let us describe the wave-function of a spin 1/2 particle as a superposition of a spin up and 
a spin down components along the direction Z of the magnetic field of the Stern-Gerlach 
device: \I/(r, t) = i/i+(r, t) \ +) + i/j-(r, t) | — ). The Pauli-Schrodinger equation describes 
the evolution of this wave-function in the presence of an external magnetic field B: 

h 2 

ihd t ^(r, t) = A^(r, t) - /xB-S^fr, t), (29) 

2m 

where S represents the Pauli matrices and is the gyromagnetic coupling constant of 
the (neutral) particle. In a Stern-Gerlach device, the field is always parallel to the Z 
direction, and its gradient is constant (B x = B y = 0, and B z = (a^ + a\z). Then, the 
evolutions of the two spin components decouple and we get that: 

h 2 

ihdtijj + /-(r, t) = -— A*0 + /_(r, t) - / + \i (a + aiz)if)+/-(r, t) (30) 

2m 

When the incoming wave-packet is initially Gaussian shaped, this equation is separable 
in Cartesian coordinates, and we find : 

h 2 

ihd t i)l h (x, t) = -—dli) x +/ _{x, t) 
ihd t ^ +/ _(y, t) = -^-d 2 ^ +/ _(y, t) 

h 2 

thd t ^ z +/ _(z, t) = -—d 2 ip z +/ _(z, t) -/ + //(a + a lZ )i; z +/ _(z, t) (31) 

The two first equations correspond to a free propagation, and in the case of an initial 
Gaussian wave packet, their solution is well known and given in many standard textbooks 
of quantum mechanics, sothat we shall not discuss how we obtain their contribution to 
the solution. The third equation is less common. We shall show, following Bohm himself 
how to solve it. Let us try to find a generalised plane wave solution of the form 
V^z-O, t) = f+/-(z, t) -expi(kz - TikH/2m). Then, / fulfills : 

ihd t f^_(z, t) = ~^(d 2 + /- 2ikd z )f*,_(z, t) -/+ M (a + a lZ )f^_{z, t) (32) 

If we could neglect the derivatives, we would have f+/_(z, t) = exp(+/ — i^(ao + a\ z)), 

but then the derivatives are not zero but yeld |^{+/ - fci^p + (^) 2 }. This term 
depends on the time only so that we can compensate it by multiplying the postulated 

2 2 

value of / by exp(i{— / + ^kt 2 — ^^t 3 }). We find so the exact solution : 

ipl/Jz, t) = exp(i(kz - Tik 2 t/2m + /- ^(a + a x z) - I + ^kt 2 - ^-^t 3 ) (33) 
' n 2m bmn 

When t = 0, this wave is plane: ip+/_(z, 0) = exp(ikz). Thanks to the linearity of the 
Schrodinger equation, the general solution is a superposition of these generalised plane 
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waves with as weight the Fourier components at time zero. The initial packet being 
Gaussian shaped, we have : 



t = o) 



| - z * 

^(773) 



Its Fourier transform is easily obtained thanks to the well-known formula : 

b 2 



[ dvexp(—av 2 + bv) = \[—exp(—) 
JveR, V a 4a 



It gives y'^fi exp(—5r 2 k 2 ) 

The wave function at time t is thus : 



^'.(z, t) = ( dkexp(i(kz -hk 2 t/2m + /- ^(a + a lZ ) ~/ +—kt 2 - ^lt 3 ) 

' J2n J n 2m 6mn 



x J*I± exp(-5r 2 k 2 ) (34) 
V V 71 

We can let go out of the integral the terms independant of k, and reorder the other ones: 

i 2 a\^ 
6mh 



r +/ -(z, t) = exp(+/-^(a + ai z)- S* 3 ) 



x 



/ dkexp(-k 2 (5r 2 + iht/2m) + ik(z - / + ^ltt t 2 )) 
J 2m 



(35) 



We get : 

r +/ -(z, t) 



,„ N exp(+/-^(a + a lZ ) - ^t 3 ) 



6mh 



x exp — 



^t 2 ) 2 

2m ' 



Using the equalities 



4(5r 2 + iht/2m) 

(org - iht/2m) 



(36) 



Sr 2 + iht/2m (5r$ + h 2 t 2 /4m 2 ) 



■ ht - 
' 2mSro - 



exp — % arccos 



Sr 



(5r 2 + 



h 2 & ■ 



and 



we get the contribution of ip+,_(z, t) to the wave-function: 



{2^Sr t y/ 2 



exp — 



(z - I + a 1 fit 2 /2mf 
A5 2 r t 
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. M (z -/ + a 1 [it 2 /2m) 2 arccos^ . fit , . ii 2 a\ ±3 . , nn . 

Notice that is obtained in exactly the same way, with the requirement that 

&o = 0,1 = 0. 

In conclusion, 

(x - v tf 



= (20F5r t )3/ 2 ex P 4 5r 2 

Z/ 2 .,hty2 arccos^ 

X 6XP 4^F' eXPH 2m5 2 r 45r 2 2 j 

(z -/+ a lf it 2 /2m) 2 .M(z -/+ a x at 2 llm) 2 arccos^- 

xexp - - — —r 1 - — — ■ expi 1 



A5 2 r t ' 2m5 2 r ASr 2 



x exp(+/ - ^ • (ao + aiz) - ^t 3 ) (38) 

with <5r 2 = 5rg-(l + t^t~j)- The interpretation of this solution which appeared originally 
in 0) is straightforward: the wave-packet of the particles of spin up (down) diffuses 
under the influence of the free part of the Schrodinger evolution (the Laplacian term) 
and is simultaneously uniformly accelerated upwards (downwards). It "feels" effectively 
a potential due to the gyromagnetic coupling which varies linearly in z, and changes of 
sign with the z-spin. 

After a time r, the particles leave the Stern-Gerlach device, and the Gaussian packets 
move freely, diffusing around centers which conserve the velocity that they possessed when 
leaving the magnets. The dependence in z is thus, up to irrelevant global phases: 



exp 



[z — I + aijirt/m + /— a 1 /ir 2 /2m) i 
4<5r 2 



2mo z r Aorf 

while the dependence in x and y is the same as before. 

In a typical single Stern-Gerlach experiment involving silver atoms, 

• the mass of the incoming atom is equal to 1,80 10~ 22 grams 

• the gyromagnetic coupling constant fi is equal to 9,27 10~ 21 gram cm 2 sec -2 Gauss -1 

• the gradient of the magnetic field along the Z axis is equal to 10 4 Gauss cm -1 

• the length of the magnet is equal to 10 cm 
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• the velocity of the incoming particle along the X axis is equal to 10 4 cm sec -1 

• the time of flight through the magnet r is equal to 10~ 3 sec 

• h is equal to 1,05 1CT 27 ergsec 

• the spreading of the incoming particle Sr is equal to 1CT 3 cm. 
Then, 

* = = 2 - 91sec " (40) 

a = ai///2m = 2, 58 10 5 cm sec -2 (41) 

P = a^/mSrl = 5, 15 K^cm -1 sec -2 (42) 

The heights of the centers of the outgoing bundles moving up and downwards after the 
time r are equal to +/— a r 2 = +/— 0,258 cm. Their Z- velocities are +/— 2a r = +/ — 
5, 15 m sec -1 . The spreading is then, taking the diffusion into account, 5r ■ vl + k 2 r 2 = 
10 -3 - \J\ + (2, 91 10 -3 ) 2 cm (approximatively 10~ 3 cm), which is very small in comparison 
to the distance between the bundles (0,515 cm), sothat these bundles are clearly separated. 
If we place a screen at a distance 1 meter after the magnet, it is reached after 10~ 2 seconds, 
and the spots are separated by a distance of order of 11 cm, while the extent of a single 
spot is still of the order of 10~ 3 cm. 



5 Appendix 2: The two particles case. 

The Pauli-Schrodinger equation allows us to describe the evolution of the wave-function 
inside the magnets of the devices: 

h 2 

ih d t tt(r L , r R , t) = - — (A L + A R )V(r L , r R , t) + (/i B L ■ S L + n B R ■ (r L , r R , t) 

2m 

(43) 

with 

A L (A R )(^ +/ _ L (r L , t) ■ Vv/-' r(tr, t) ■ | + / - ) L (8) | +' / -' ) R ) = 
(A L ^ +/ _ L (r L , t)) ■ V'+v-' *(r fl , t) ■ | + / - ) L (8) | +' / -' ) R 
(^ +/ _ L (r L , i) • (A fl V + 7-' 0) • I + / " )l ® I +7 >r) ( 44 ) 

where A £ / fi is the Laplacian operator in the Left/Right coordinates, B L / R is the 
magnetic field in the Left/Right regions, while the components of Sl/r are the Sigma 
matrices of Pauli. When the magnetic fields ares parallel to the Z axes, only the third 
Sigma matrices appear, which are defined by: 

S 3 (J?) ^ +/ _ L (r L , t) ■ VV/-' r(t r , t) ■ | + / - ) L (g) | +' / -' ) R = 
+ / " (+7-V+/- L(r L , 0) ' V'+V-' fl(r fl , *) • I + / " >l ® I +' / >r («) 
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Then, it is easy to check that the Pauli-Schrodinger equation is separable into the six spa- 
tial coordinates and that the general solution, when the wave-packet is initially Gaussian, 
can be expressed in terms of the solutions associated to the single particle case that we 
obtained in the previous appendix. 

Let us now derive the conservation equation associated to this evolution law. We can 
write it in the form: 

h 2 

ihd t V(r, t) = -— A*(r, t) + /iB • £*(r, t) (46) 

where the Laplacian is a 6-dimensional Laplacian, the magnetic field a 6-dimensional field. 
It is worth noting that £ is a self-adjoint operator on the spin space (isomorph to C 4 ). 
The adjoint (transposate-conjugate) of the previous equation is: 

h 2 

- ih d t * f (r, t) = A^(r, t) + * f (r, t) p B • S (47) 

2m 

where we used the fact that £ = £L Let us now multiply the equation of evolution by 
\frf, its adjoint equation by \I/ and take their difference, we get: 

h 2 

ihd t * f (r, t) ■ tt(r, t) = - — * f (r, t) ■ (Atf(r, t)) - (A^ f (r, t)) ■ (r, t) (48) 



We can rewrite the second member of the last equation as the 6 dimensional divergency 
of a 6 dimensional flux: 

#t(r, t) ■ (A*(r, t)) - (A&(r, *)) • #(r, t) = 

div(^(r, t) ■ grader, t) - grad* f (r, t) ■ *(r, t)) (49) 

Now, ^L(r, t) • grad\E'(r, t) is the complex conjugate of grad^L(r, i) • \&(r, t) so that we 
obtain the equation of conservation d t p = div(J) where: 

p = ¥(r, t) ■ *(r, t) (50) 

or, in another notation: 

p(r L , r R , t) = (tf(r L , r^, t) |tf(r L , r fl , t) ) 

and 

J = (h/m)S(¥(r, t)grad^(r, *)) (51) 
where $$(z) the imaginary part of z. 

The three and four particles cases can be treated in a similar way. 
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